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The energy structure of the V15 is studied using exact diagonalization methods, and the 
adiabatic changes of the magnetization in the system with the sweeping field are investigated. 
We confirm that the Dzyaloshinskii-Moriya interaction leads to an energy gap which allows the 
adiabatic change of the magnetization which has been found experimentally by Chiorescu et 
al., and also predict novel dynamics of the magnetization due to this energy level structure. 
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§1. Introduction 

Interesting adiabatic change of magnetization in V15 was discovered by Chiorescu, et al., 1 ) because 
of which the energy level structure of V15 has attracted much attention. Detailed information of the 
energy levels is necessary in this case, since time dependence of magnetization is almost completely 
governed by pure quantum dynamics. Besides, we can also expect other interesting phenomena 
related to quantum dynamics in this material, due to interesting energy level structure that could 
emerge from new models. 

The time-reversal symmetry and the half-odd value of the magnetization imply that the eigenstate 
at zero field must be at least doubly degenerate and therefore the quantum tunneling may be 
suppressed. However, this idea contradicts the observed adiabatic change of the magnetization in 
V15 in so far as we consider only a single mode of magnetization. If the geometric frustration of 
interactions in V15 is taken into account, one obtains a doubly degenerate avoided level crossing 
structure. 2 ' 3 ^ That is, two avoided level crossing structures exist and they cross each other at H = 0. 
The adiabatic change of the magnetization occurs in each avoided level crossing structure. 2 ) This 
degenerate structure has brought forward an interesting question: how does the magnetization 
change from 1/2 to 3/2? It has been shown, using a three spin model 3 ) that one of the states with 
M = 1/2 evolves to a state of M = 3/2 while the other state remains in the same M = 1/2 state. 
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Recently more realistic calculations have been done on the model which includes all the 15 spins 



The V15 molecule has a C3 rotational symmetry and lacks a symmetry plane passing vertically 
through the middle of the molecule. That is, V15 molecule does not have the symmetry between 
the right-hand rotation and left-hand rotation. This implies that the DM interaction in the xy 
plane of the molecule can be non- vanishing. The DM interaction causes mixing of the two pairs of 
low energy states of S = 1/2 to form two sets of avoided level crossing states as a function of the 
external field H. Thus the system has an energy gap at H = and shows an adiabatic reversal of 
the magnetization. Konstantinidis and Coffy 5 ) studied the energy structure of the system with the 
DM interaction by a perturbational method and obtained the dependence of magnetization on the 
strength of DM interaction. 

We extend their model by exploiting the C3 symmetry and calculate the energy structure by 
exact diagonalization method. 6 ) We employ the full 15-spin Hamiltonian with DM interactions, 



of Vi5(Fig.l). 4 ' 5 ) Experimentally there is evidence of a small splitting between the lowest doublets 
of V15. It is believed that Dzyaloshinskii-Moriya (DM) interaction is the cause of this gap. 2,3,5 ) 
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Fig. 1. Structure of V15: the bold line denotes J, the thin line J2, and the dashed line Ji 
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§2. Method of Computation 

We have treated the DM interaction as a perturbation over the eight low-lying states of the 
Ho- The low- lying energy levels of Ho are very closely spaced with many degeneracies. To obtain 
the eigenvalues and eigenvectors with high accuracy and to gain information about the nature of 
the states we have used C3 symmetry to divide the space into the doubly degenerate E and the 
nondegenerate A subspaces. We set up the Hamiltonian matrix of the Heisenberg Hamiltonian Hq 
in the symmetrized basis set in the chosen subspace. We obtain the eigenvalues and eigenvectors 
of the lowest eight energy levels using Davidson's algorithm. 6 ) In the next step, we set-up the 
perturbation matrix Hbm in the basis of the eight low- lying eigenstates of Ho- This is achieved 
by setting up the full matrix of Hbm in the full unsymmetrized basis and then transforming the 
matrix into the basis of the eight low-lying eigenstates. To check the accuracy of the method, we 
have also diagonalized the total Hamiltonian calculated in the full unsymmetrized basis (where DM 
interaction is treated in a non-perturbative manner) with D x = D y = D z =40 K and compared with 
the results from the above approximate method. We found not only the eigenvalues but also the 
eigenvectors of the lowest few states obtained by both methods to agree with each other very well 
for the lowest few eigenstates that could be obtained from the full Hamiltonian matrix. Therefore, 
we use the former method to calculate the quantities discussed in the following. 

Using various sets of the parameters (J, Ji, J2, Dij), we have fitted the experimentally known 
energy level structure. In order to fit the energy gap between the ground state doublet and the 
first excited quartet in Ho, we adopt the parameters 4 ) 

J = 800K, Ji = 54.4K, J 2 = 160K. (2.1) 

The configuration of the DM vectors is not known yet. At this moment we make use of the 
fact that the zero field energy gap of the avoided level crossing is about 80mK, 7 ) and the width of 
the magnetic jump from M z = 1/2 to M z = 3/2 is « 0.75K. 8 ) Although we have only two values 
to be fitted, there are many bonds in the molecule and therefore many DM parameters. To reduce 
the set of DM parameters, we follow Konstantinidis and Coffy and assume that the vector 
is nonzero only for the pair of spins between which the Heisenberg exchange J is strongest. We 
further assume that the DM vectors of the upper hexagon and lower hexagons to be the same to 
minimize the number of free parameters in the model. 



3 




4 




D 



Fig. 2. In all the figures Gap A (unit of Kelvin) is plotted on left Y axis and Width (unit of Tesla) on right Y 
axis. In Fig. (a) D x = D y — D z . In Figs, (b), (c) and (d), we set D y = D z — 0, D x — D z = and D x = D y — 0, 
respectively. D x =10 K in Figs, (e) and (f). D y =W K in Figs, (g) and (h). D z — 10 K in Figs, (i) and (j). In Figs, 
(e) to (j), the component of DM vector, other than the one taken as :r-axis and the one fixed at 10 K, is set to 
K. D x = D y ; D z = in Fig. (k). D x = D z ; D y = in Fig. (1). D y = D z ; D x =Qm Fig. (m). 

We have explored the parameter space of Dij quite extensively and show in Fig. 2 the variation 
of the energy gap between the two doublets and the width of the transition at the level crossing 
between M=l/2 and M=3/2. There, (D x , D y , D z ) means the DM vector on the bond of J, say 
Di2- The DM vectors at other two bonds of a hexagon are given by a rotation of 120°: 

£>f 4 = \ (~Dx ~ VSD y ) , D y u = \ (VSD X - D y ) , D* 34 = D z 

L»f 6 = \ (-D x + VSD y ) , D y m = \ (-VSD X - D y ) , D% % = D Z 

We note that the gap depends strongly on D z as we see in Fig.2 (d), (f), (h), (1) and (m), and 
the width strongly on the components D x and D y as we see in Fig.2 (b), (c), (e), (i), (j) and (k). 
Strong dependence of width on D x and D y can be understood from the fact that they contain 
terms which change the M z value of the state whereas D z part conserves the M z symmetry. We 
found the set (D x , D y , D z ) = (40,40,40) is a good candidate to reproduce the experimental gap 
and width. 
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We find the energy structure is similar to that of Model V in our previous study. 3 ) In Fig.3, the 
eight low-lying energy levels are depicted as a function of the field. Around H = the low lying 
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Fig. 3. Energy levels obtained for Vib (a) for the set of levels S=l/2 and S=3/2. In (b) the avoided crossing between 
the M=l/2 and M=3/2 levels is shown on a finer scale. 

four states forms two sets of avoided level crossings. They are almost degenerate at small fields. 
However at the critical field Hq where the ground state magnetization changes from 1/2 to 3/2, 
one of the degenerate states adiabatically changes to a state of M z ^3/2, while the other does not 
change the magnetization. This structure is a consequence of the double-degeneracy of the S = 1/2 
ground states and causes interesting adiabatic process as discussed in the next section. 

We find that a model with three spins interacting antiferromagnetically has the same symmetry 
of V15, and the model gives a good approximation for the present model if we choose J = 2.44K 
and D y = 0.06 J. 

§3. Adiabatic change of the magnetization 

The magnetization process of the present model shows a peculiar adiabatic magnetization process 
in a sweeping field. Because the two avoided level crossing structures are very close to each other 
near the zero field, we assume that the two low energy states are equally populated. When the 
field is swept, one of them changes to a state of M = 3/2 and the other keeps M = 1/2. Thus, in 
the adiabatic motion, the magnetization changes as 



observed at low temperature where the thermalization process is slow, although the magnetization 




(3.1) 



instead of |. 



We, therefore, expect that the magnetization changes to 1. This change will be 
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1 finally relaxes to | due to the thermal effect. If we sweep back the field, the magnetization 
process will show a loop structure due to the thermalization effect. The process in a dissipative 
environments will be studied elsewhere. 9 ) For very slow sweeping, the thermalization causes the 
system to realize the equilibrium magnetization process and the magnetization directly changes to 
M = 3/2. 

§4. Hidden conservation law 

Here we would like point out another characteristic feature of the energy level structure. As we 
have pointed out in a simple model, 3 ) there exists a kind of selection rule. If we adiabatically change 
a state starting from the ground state at large negative field where the magnetization is almost 
—3/2, the state stays in the ground state until H = 0. On the way the magnetization adiabatically 
changes from —3/2 to —1/2. At H = the state crosses with another state. Near H = 0, both 
states almost are degenerate and it is difficult to see the difference. The magnetization changes 
from —1/2 to 1/2 around H = 0. When we increase the field further, in the adiabatic change of 
the state, the magnetization does not change at the critical field Hq where the magnetization of 
the ground state changes from 1/2 to 3/2. That is, the magnetization adiabatically reached is 1/2 
for H = +oo, although the ground state magnetization is 3/2 for H = +oo. 

Here we adopt the following method to study the quantum dynamics of the model. We construct 
the total matrix 7i as mentioned previously in the basis of lowest eight states but the field is now 
time dependent. We choose the ground state at large negative field of —6 Tesla as the initial state 
|Y>o >■ Then, the field is swept at a speed v for a time t up to a large positive field of 6 Tesla. We 
have studied the time evolution of the system by solving the time-dependent Schrodinger equation 

ih^=HW. (4.1) 

The state is time-evolved by the equation 

1>{t + At) = e -<W(*+^)At/» ^) . (4.2) 

The evolution is carried out by explicit diagonalization of the Hamiltonian matrix H(t + 
and using the resulting eigenvalues and eigenvectors to evaluate the matrix of the time evolution 
operator e ~*' w ( i +~) At / 7i . Here we take At = 1CP 5 , which is small enough because the maximum 
eigen value of the energy is E max = —4.598 (in unit of J) and E mSLX At <C O(l). We sweep the 
magnetic field with the velocity v = 2.0 x 1CT 6 , which is slow enough for the adiabatic change, 
because the energy gap at H = is 80 mK and thus the Landau-Zener-Stiickelberg 10 ) transition 
probability is estimated 

p = l-exp(- 7r( ^ ) ) = 0.007. (4.3) 

In Fig. 4, the adiabatic change of the magnetization as a function of H is depicted for the 15 spin 
model. This property is also found in the corresponding triangle model. 
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Fig. 4. Adiabatic change of the magnetization starting from the ground state with M z — —3/2. 



§5. Summary and Discussion 

We have investigated the energy level structure of V15 using exact diagonalization method. We 
find various peculiar properties of the energy level structure which are also expected from earlier 
studies of an effective three spin model. 3 ) First of all, we have shown that the DM exchange interac- 
tion leads to an energy gap which allows for adiabatic change of the magnetization. ^ Furthermore, 
the two-fold degeneracy of the avoided level crossing structure leads to an unusual energy level 
ordering in which the ground state magnetization changes from 1/2 to 3/2. 

Because the adiabatic change is explicitly observed experimentally, we expect that the novel dy- 
namics of the magnetization due to this energy level structure can also be discovered in experiments. 
In this way we hope we can establish the role of real time quantum dynamics in magnetic systems, 
although we need a rather fast sweeping of the magnetic field in order to avoid the dissipative 
effects. 11 ) Studies of the dissipation effect will be reported elsewhere. 9 ) 
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